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Abstract. Let K be the product 0(m) X 0(n 2 ) x ••• x 0(n r ) of orthogonal groups. 
Let V = <8>i =1 C ni , the r-fold tensor product of defining representations of each orthogonal 
factor. We compute a stable formula for the dimension of the if-invariant algebra of degree d 
1 homogeneous polynomial functions on V. To accomplish this, we compute a formula for the 

number of matchings which commute with a fixed permutation. Finally, we provide formulas 
for the invariants and describe a bijection between a basis for the space of invariants and 
the isomorphism classes of certain r-regular graphs on d vertices, as well as a method of 
associating each invariant to other combinatorial settings such as phylogenetic trees. 

1. Introduction 

Let r G Z + , d G Z> , and let n = (ni,n 2 , . . . , n r ) be a sequence of positive integers. Let 
V(n) = C ni <g> C na <g> ■ ■ ■ ® O and G(n) = GL ni (C) x GL n2 (C) x • • • x GL Ur (C), where 
GL ni (C) denotes the general linear group of degree over the complex numbers. Finally, we 
denote by K(n) the product O(rii) x 0{n 2 ) x • • • x 0{n r ), where 0(rii) denotes the complex 
orthogonal group of degree defined by 0{rii) = {56 GL m (C) \ g T g = I}. The group 
K(n) acts on V(n) by 

(^1, 92, ■ ■ ■ , 9r) ■ {v v , v 2 ,..., v r ) = (gxvx) <8) (g 2 v 2 ) ® • • • <g> (g r v r ) 

for j.GOfn.j.D.GC"'. 

Let V(V(n)) denote the algebra of complex valued polynomial functions on V(n). The 
group K(n) acts on V(V(n)) in the standard way: for k G K(n), v G V"(n), and / G 
V(V(n)), we have k ■ f(y) = f(k T v). The i^(n)-invariant subspace of V(V), which we 
denote by V(V(n)) K( - n \ is known to be finitely generated [|fiii90 |; however, a description of 



these invariants is incomplete outside of a few particular values of n. 

We have the standard gradation V(V(n)) = 0^ o V d (V(n)), where V d (V(n)) is the 
subspace of degree d homogeneous polynomials on V(n). Moreover, the i^(n)-invariant 
subalgebra inherits this gradation. That is, V d {V{n)) K{ - n) = V d {V(n)) n V(V(n)) K{n \ 

While dimV d {V(n)) K ^ is unknown in general, we show in section |2| that the dimensions 
stabilize as the n^'s are sufficiently large. It is in this setting we seek a formula for these 
dimensions, as well as a description of the invariants themselves. We first note that if d is 
odd, we immediately obtain 

dimP d (V(n)) A ' (n) =0 

and so we assume d = 2m where m is a positive integer. 

We call A a partition of size d, denoted A h d, if A = (Ai > A2 > • • • > \e > 0) is a weakly 
decreasing and finite sequence of positive integers such that Yli.=i = d. If A has b\ ones, 
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&2 twos, bs threes, etc, we say that A has shape (l bl 2 b2 3 &3 • • • )• We define 



l^ 2 b ^ ■■■b 1 \b 2 %\--- 



We let £(X) denote the length of the partition. Finally, we say A = (Ai, A 2 , . . . , A^) is even if 
A; e 2Z for all i. 

Let a be an element of S m , the symmetric group on m letters. When written in disjoint 
cycle notation, the lengths of the cycles of a form an integer partition of m, and we refer to 
the shape of a as above. We call a permutation r a matching if it contains only cycles of 
length two; that is, r is a fixed point free involution. 

In section |[ we show that 

Theorem 1. Let d be a positive, even integer and let K(n) = 0(ni) x 0(ri2) x • ■ ■ x 0(n r ) 
and V(n) = C ni ® C n2 <S> ■ • ■ <%> C n '', where ni, n 2 , . . . ,n r > d. Then dimension of the space 
of K(n) -invariant degree d homogeneous polynomial functions on V(n) is 



dimP d (K(n))^ (n) = 



Xhd 



N(\y 



where N(\) is the number of matching s that commute with a permutation of shape A. 

We determine the number of matchings which commute with a given permutation in 
section 0. In particular, we have 



Theorem 2. Given a permutation a with shape A 
that commute with a is given by 



'l bl 2 b2 ■ ■ -t bt ), the number of matchings 
N{\) = iV((l bl )) • iV((2 62 )) N((t bt )) 



where 



N((a b )) 



b\a b / 2 
2 ft /2(|)! 



E 



i=l, i odd i\(tz±')\2(b-i)/2 



6!a( 6 - i )/ 2 



i even ji^ 2=1 ynfb-i) /2 



if a odd and b odd 
if a odd and b even 

if a even and b odd 

if a even and b even 



The notion of matchings has surprisingly deep connections to a variety of areas of research. 
Classical Schur-Weyl duality describes a relationship between irreducible finite dimensional 
representations of the general linear group and the symmetric group. The natural actions 
of each of these groups on the tensor space <§Z) k (C n ) centralize each other, resulting in the 
multiplicity free decomposition 



^k C n ^ F X q w > 



Xhk 
e(X)<n 



where F x and W x are the irreducible representations of GL n and S n , respectively, associated 
to the partition A. In 1937, Richard Brauer defined an algebra which plays the role of the 
symmetric group in a similar statement on the representation theory of the orthogonal group 
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Pra37|| . For a partition A = (Ai, A2, • • ■)> ^ -V = (A' 1; A' 2 , . . .) denote the conjugate partition 
(thus X'j is the number of boxes in the jth column of the Young diagram of A). We have 

[k/2] 

® k c n ^0 v x ® u x 

i=0 X\-(k-2i) 
X^+X'^n 

where V x and U x denote the respective irreducible representations of the Brauer algebra 
£>fc(n) and 0{n). The dimension of the algebra Bk{n) is {2k — l)(2k — 3) ■ ■ ■ 3 ■ • • 1, the 
number of matchings on 2k elements. These basis elements of the algebra are frequently 
depicted as graphs; the diagram below represents a possible basis element of the algebra 
B 6 (n): 



The number of matchings on a set of elements has applications outside of representation 
theory as well. Numerous fields in biology use graphs called phylogenetic trees to illustrate 
inferred evolutionary relationships. These graphs are full rooted binary trees in which each 
node, with the exception of the leaves, has exactly two children. Each node could, for 
example, represent a particular species. A parent node (that is, one that is not a leaf) would 
represent the most recent genetic ancestor of the two child species. Diaconis and Holmes 
PH98| have described a bijection between matchings on 2n elements and phylogenetic trees 
with n + 1 leaves, which we recount here. 

We first describe the procedure for building a tree with n+1 labelled leaves from a matching 
on the set {1, 2, ... , 2n}, which we illustrate with the example r = (1 4)(2 3)(5 8)(6 7). The 
resulting tree will have leaves labelled with the set {1, 2, . . . , n + 1}; the ancestors of the tree 
will be labelled {n + 2, n + 3, . . . , 2n}. First, look for any pairs in the matching that contain 
only numbers at most n + In our specific example, we have two such cycles, (1 4) and 
(2 3). Each such cycle defines a sibling pair of leaves. Choose the pair with the smallest 
child, in this case, (14). Label the parent of this pair with the smallest ancestor, 6. The 
other sibling pair will be labelled with the next available ancestor, 7. Since 6 and 7 appear 
in the same cycle of the matching, they must also be siblings in the tree. Their parent will 
be the last available ancestor, 8. Finally, 8 must be paired with 5, which will form its own 
leaf. The resulting tree is shown on the left below. The corresponding tree with labelled 
leaves is on the right. 




1 4 2 3 1 4 2 3 
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Note that there is more than one possible rule to create such a correspondence. 

We now sketch a method of finding a matching when given a tree. We start with a binary 
rooted tree with n + 1 labelled leaves. We look for the pair of siblings with the smallest child, 
and label the parent of these children with n + 2. Repeat this process until all nodes (except 
the root) have been labelled. The matching defined by this tree is formed by pairing siblings 
into cycles. A particular example is shown below. The tree on the left is a binary tree with 
7 leaves. The center diagram is obtained by labeling the nodes of the tree as described, and 
the corresponding permutation is on the right. 




(1 8) (2 7) (3 10)(4 11)(5 6)(9 12) 



5 6 5 6 



In section f|, we show that given an r-tuple of matchings (r 1; r 2 , . . . , r r ), each on 2m letters, 
we can describe an invariant in V 2m (V(n)) K ( n \ The symmetric group S^m acts on such a 
tuple by simultaneous conjugation; that is, for g G S2 m , we define 

g • (n, r 2 , . . . , T r ) = (gng- 1 , gr 2 g~ l , gr r g~ l ) 

We show that the orbits of this action are associated to a generator of the algebra V 2m (V(\\)) K(ji \ 
and are in bijection to the isomorphism classes of certain graphs. 

We can therefore consider an action of S2 m on r-tuples of phylogenetic trees with m + 1 
leaves. Using the method described above, we write an r-tuple of trees as an r-tuple of 
matchings, apply the action of simultaneous conjugation by an element of 5*2 m , and then 
draw the list of trees associated to the resulting r-tuple of matchings. Hence, we define an 
action of S2 m on a forest of r phylogenetic trees, each with m + 1 leaves. 

To illustrate, suppose we choose a tuple of three trees, each with four leaves. Labeling the 
roots of the trees 1 through 3 and adjoining these roots to a common vertex creates a forest: 




13 2 4 



By following the process outlined earlier, we associate this forest to the 3-tuple of matchings 
(Ti,T 2 ,T 3 ) = ((1 3) (2 5) (4 6), (1 3)(2 4)(5 6), (1 6) (2 4)(3 5)) 
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We choose a permutation g = (1 3 5) (2 4) (6) to act on this tuple: 

(gng- 1 ^^- 1 ^^- 1 ) = ((1 4) (2 6) (3 5), (1 6)(2 4) (3 5), (1 5)(2 4)(3 6)) 
and draw the forest associated to the result: 




In [ HWW11 |, a bijection is described between the invariants defined on the tensor space 
®[ =1 C™ 1 under the action of a product of unitary groups and isomorphism classes of finite 
coverings of connected simple graphs. In section f|, we seek a similar graphical interpretation 
of the invariants on the same space under the action of orthogonal groups. Specifically, we 
produce formulas for the invariants as well as a bijection between a basis for the space of 
invariants and isomorphism classes of edge-colored r-regular graphs on d vertices. 



2. Proof of Theorem [1] 

Let m be a positive integer and let d = 2m. For each partition A of d with at most n 
parts, there exists an irreducible representation of GL n (C) with highest weight indexed by 
A, which we denote by F£ following [ GrW09 |. If A(i) h d for 1 < i < r, then 

is an irreducible representation of G(n) which embeds in V d (V(n)) with multiplicity denoted 
by 9\m\(2)...\(r). That is, 

V d (V(n)) = SaU)a( 2 )..^< (1) ® ® " • • ® ^ 



We can now write the f^(n)-invariants as 



[7> d (V(n))]* (n) - , A(1)A(2 ,.. AM (F r f ® F„f > ® . . . ® Ff ' 



A( 4 )hd 
^(A(0)<ni 



A(D \ ° (ni) ^ / r A( 2 ) \° (n2) ^ ^ / cAW \ ° ( "' r) 



^(A(^))<ni 



The Cartan-Helgason Theorem ||GW09|| tells us that if F x is an irreducible representation 
of GL n (C), then dim (F x )° (n) is at most one. In particular, we have 

dimfi?) ^/ 1 ifAiseven 

I otherwise 

and so 

dim[V d (V(n))] K{n) = 0ACDAW...AC) C 1 ) 

A(i),...,AW 

A(0|-2m 
AM even 

The dimensions will stabilize when the sum is taken over all even partitions of 2m; that is, 
when ni, n 2i ■ ■ ■ j n>r > 2m. 

We next recall Schur-Weyl duality, which will allow us to find another interpretation of 
Equation |l| While GL n acts on the space cg> m C n by simultaneous matrix multiplication, 
the symmetric group acts on the same space by permuting tensor factors. That is, given 
x G GL n and a permutation a G S m , we have 

x(vi ® v 2 ® ■ ■ ■ <g) v m ) = xv i g) xv 2 <g) • • • <g) xw m 

and 

Cr(t>i ® -U 2 <8) • • • <8) -U m ) = D CT (1) <g> W CT ( 2 ) <g> ■ ■ ■ <8> W<T(m) 

for vi,V2, ■ ■ ■ ,v m G C n . The subalgebras of End(® m C n ) generated by each of these actions 
are full mutual commutants. Hence, we obtain a multiplicity free decomposition 



where f/^ is an irreducible representation of S'm indexed by the partition A. 

Let U x denote the irreducible representation of the symmetric group SVm indexed by the 
partition A. We have 

dim^ 1 ) ® f/ A ( 2 ) <g) • • • ® f/ A «)^ = ^ (1)A(2) ... AW (2) 

where each A(i) h 2m. These are in fact the same multiplicities that appear in the de- 
composition of representations of the general linear group | |HW09| , Prop. 3]. The numbers 
fi'A(i)A(2) --A(r) are known as the Kronecker coefficients, and are known to be hard to compute. 

Fix r = (1 2) (3 4) • • • (2m — 1 2m) in S^m- Denote by H m the centralizeij] of To in S2 m - 
We note that (S 2 m,H m ) is a Gelfand pair [HSS92|| , and so dim(U x ) Hm is at most one. In 



fact, we have 

1 if A is even 
otherwise 



dim{U x ) Hm 
Thus, 

dim[[/ A(1) (g) U m <g> ■ ■ ■ <g) u x W] HmXHmX ""* Hm 



1 if is even for all i G {1, 2, . . . , r} 
otherwise 

(3) 



1 -ff m is isomorphic to the hyperoctahedral group, the group of symmetries of a hypercube of dimension 

rn. 
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As an immediate consequence of Peter- Weyl decomposition, we have 

C{S 2m ] = Q)U x ®U x 

A 

We can now describe the group algebra of a product of r copies of S 2m as follows: 
C[S 2m xS 2m x---x S 2m ] C[S 2m ] ® C[S 2m ] ® • • • <g> C[S 2m ] 



A(!),...,AM 

Combining this fact with Equations and ^J, we obtain 
dimC[S 2m x S 2m x • • • x 5 2m ]5 2ro x^ m x...xif m 



U x <g> f/ A j <g> f f/ A ® £/ A j <g> ■ • • ® f £/ A ® f/ A 

(u x(1) ®u x<1) )®---® (u x<r) ® u x(r) ) 



dim ^ 



f/ A(1) ® • • • ® U xM ) S2m ® (U XW ® • • • ® f/ AWX 



A(i),...,AM 
AWh2m 

= ^AWaW-AM 
AW,..., AM 
A( i )h2m 

aW even 
So far, we have shown that 

dim [P d (V(n))] K(n) = dimC[S 2m x S 2m x ■ ■ • x S3m ] s *»* H <»*"* H >> 
Define the set 

hm = {r <E S 2m : r 2 = id, r(i) ^ i for all % < n} 

That is, I 2m is the set of matchings on 2m letters. As an S^-set, we have I 2m = S 2m /H m . 
We will denote the product of r copies of I 2m by T^ m . 

Recall that given a group G and H C G, we have dimC[G]^ = dimC[G/H}. Thus, by 
the previous results 



^ 9\W\W...\(.r) = dimCfS 1 ^ x S 2m x • • • x S>. 



, ^ H m X ■ ■ ■ X //rr 

2m 



A(!),...,AM 
A«|-2m 

AM even 

= dimC[S 2m /# m x • • • x S 2m /H m ] s *™ 
= dimC[TU S2m 

Given a group G and a set A, it is well known that the dimension of the space of invariants 
C[A] is equal to the number of orbits of the action of G on A. By Burnside's Lemma, the 
number of such orbits is the average number of x G A fixed by g G G; that is, 



dimC[A] G = tLJ2\ X 



geG 



where \X\ 9 denotes the cardinality of the set of points in X fixed by g. In our setting, where 
G = S 2m and X = /J™, we have 



\02m.\ 



where S 2m acts on /J m by simultaneous conjugation. That is, given (n, r 2 , . . . , r r ) G / 2 r m , 
# ^ SWru we define the action 



g.(r 1 ,T 2 ,...,T r ) = (gng ■■ 1 — ■■ 1 



gr 2 g gr r g 



We easily see that 



\ir 



2ml — \\ 1 2m 



and so, given (7 G S^m; h remains only to find a formula for the number of matchings t G -^2m 
such that grg^ 1 = r. It is easily seen that this number is the same for two permutations with 
the same cycle type, as we can simply relabel the entries of each cycle. Thus, if g, h G S 2m 
have the same cycle type, we have \I 2m \ 9 = \hm\ h - 

Recall that two elements in S2 m are conjugate if they have the same cycle type. Denote 
by S 2m the set of conjugacy classes in S 2m , indexed by integer partitions /1 of 2m. Hence we 
can define a class function N : S2m — > N by setting N(X) equal to the number of matchings 
t G /2m that commute with a permutation g G S2m with cycle type fx. Finally, we have 
shown that for V(n) = C ni ® C n2 ® ■ ■ • <g> C n '', fsT(n) = O(ni) x 0(n 2 ) x ■ • • x 0(n r ) and 
rij > 2m for 1 < i < r, we have 

dim [V 2m (V(n))] K(n) = dimCfc]*- = E l^'ml 9 ~E iV ( A ) r 

The formula for N(X) is presented in Theorem |[ which we prove in the next section. 



3. Proof of Theorem g 

We now determine the formula for N(X), the number of matchings which commute with a 
permutation with cycle type A. Our strategy will be to focus initially on permutations with 
shape (m 2 ). We then consider "brick" permutations with shape {a b ), where ab = 2m, and 
finally we generalize to all permutations. 

3.1. The Case A = (m 2 ). We begin with: 

Lemma 3. Let g = (oti a 2 ■ ■ ■ a m )(/3i (3 2 ■ • ■ (3 m ) be a product of two cycles of equal length 
m. If a is a matching that commutes with g, then either a = g m ^ 2 or o has form 

a = («! 0j)(a 2 (3 j+1 )(a 3 (3 j+2 ) ■ ■ ■ (a m 

where < % < m. Hence, if m is even there are exactly m + 1 matchings that commute with 
g, and m such matchings if m is odd. 

Proof. Suppose first that a permutes the elements within each cycle of g. Let a(c(k) = ct m 
for some k with 1 < k < m; that is, assume the transposition (a& a m ) appears in a. Then 

«fc+i = g{ak) = °~go-(a k ) = ag(a m ) = a(a x ) 

8 



and so (cti a k+ i) is a transposition in a. Similarly, 

= g(a k +i) = crgcr(a k+1 ) = ag{a 1 ) = a(a 2 ) 

and so (a 2 is a transposition in a. Continuing in this way, we find a contains the 

transpositions («i a k+1 ), (a 2 a k+2 ), (a 3 a k+3 ), ...,(a k a rn ). Now 

a(a m )a k = g{a k _i) = aga{a k _ 1 ) = ag(a 2k -i) = o{a 2k ) 

and hence m = 2k. Then we have only one possible choice for k, k = m/2. Thus 

a = (ai arn +1 )(a 2 a™ +2 )(a 3 a™ +3 ) ■ ■ ■ {arn « m ) = ( ai a 2 ■ ■ ■ a m ) ml2 

Thus if a permutes the elements within the cycles of g, we must have a = g m ^ 2 . Note such 
a a only exists if m is even. 

Now suppose a interchanges elements between the two cycles of g. Assume o[a\) = fa for 
some j with 1 < j < m. Then (aii fa) is a transposition appearing in a. Since ago = g, we 
have 

«2 = gM = crga(ai) = ag(fa) = cr(f3 j+1 ) 
and so (a 2 fa+i) appears in a. Similarly, 

«3 = g(a 2 ) = aga(a 2 ) = ag(fa +1 ) = a(fa +2 ) 
and so (a 3 fa+ 2 ) appears in a. Continuing in this way, we have 

a = («i fa)(a 2 fa+i)(a 3 fa +2 ) ■ • ■ (a m fa-i) 

□ 

The matchings that will commute with a fixed permutation with two cycles of equal length 
can easily be interpreted as diagrams. With g as in the above proof, we draw two rows of 
m nodes. We label the nodes along the top row with aii, a 2 , . . . , a m , and along the bottom 
row with fa, fa, ... , f3 m . Drawing edges so that each node is connected to exactly one other 
node defines a matching. When we consider only the diagrams corresponding to matchings 
which commute with g, an immediate pattern emerges. If m is even, the first such diagram 
is obtained by drawing an edge from a± to arn +1 , an edge from a 2 to ar± +2 , and so on. The 
nodes labelled fa are matched in an identical way, resulting in a diagram which represents 
the matching g m ^ 2 . The remaining m diagrams (for odd or even values of m) are obtained 
by drawing edges from «j to fa for each i, and then cyclically permuting the second row of 
unlabeled nodes. The five matchings that commute with g = (ai a 2 a 3 a 4 )(/3i f3 2 fa fa), for 
instance, are: 




a l &2 ®3 ^4 <%1 &2 a 3 ft 4 




fa fa fa fa fa fa fa fa 
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3.2. The Case A = (a b ). Suppose now that our permutation g has b cycles of length a, 
where ah = 2m. Recall that we can conjugate g by a permutation a by applying a to each 
symbol of each cycle of g. That is, if g = (ai ol% • • • a a )(/3i fa • • ■ Pa) • • • , we have 

aga' l = (a(a 1 ) a(a 2 ) ■■■ a{a a )){a{fa) a((3 2 ) ■■■ a((3 a )) ■ ■ ■ 

It follows that if a commutes with g, conjugation by a sends one cycle of g to another cycle 
of g. Thus we can take pairs of cycles of g and look for the matchings on 2a numbers that 
commute with the product of these pairs following Lemma |3|. If we have an odd number of 
cycles, the unpaired cycle will commute with a power of itself, that power being a/2. The 
product of these will be a matching which commutes with g itself. 

To illustrate, we'll quickly compute a matching which commutes with the permutation 
g — (1 2 3 4) (5 6 7 8) (9 10 11 12). Suppose we first pair the first and third cycle, and look 
for matchings which commute with (1 2 3 4) (9 10 11 12). By Lemma [| we have five to 
choose from, such as (1 10) (2 11) (3 12) (4 9). We still have one cycle of g left, and there is 
only one matching to commute with it: (5 7) (6 8). Hence, a matching that commutes with 
g is (1 10)(2 11)(3 12)(4 9) (5 7) (6 8). 

Once again, we turn to diagrams to simplify. After fixing a cycle expression for g, we draw 
a diagram with b rows of a nodes, and label the nodes of the ith row with the entries of the 
ith cycle of g, as before. Again, we can define a matching by drawing an edge between pairs 
of nodes. The following diagrams illustrate the case where g has cycle type A = (4 5 ). The 
diagram on the left displays a matching, but not a matching that will commute with the 
permutation g, as it fails to map one cycle of g to another. The center diagram meets this 
requirement, but does not satisfy Lemma Finally, the diagram on the right corresponds 
to a matching that will commute with g. For simplicity, we have suppressed labeling of the 
nodes. 




We have now established a convenient way of establishing which matchings commute with 
a permutation of shape (a fe ), it is simply a matter of counting them. To do this, we will 
compress the diagrams which correspond to an eligible matching as follows: for each cycle 
of the permutation, draw a single node. By Proposition [|, there are a matchings that will 
commute with two distinct cycles. Assign a color to each of these choices, and color the edge 
between the corresponding nodes accordingly. If matching contains the a/2-th power of a 
cycle, the corresponding node is left isolated. Thus N((a b )) is obtained by counting graphs 
of the form 

10 



We consider three sub cases, depending on the parity of a and b (recall we cannot allow 
both a and b to be odd, as we require ab to be even): 



3.2.1. a odd, b even. If a is odd, any matching that commutes with g cannot contain a power 
of a cycle of g, so we have no isolated nodes. The number of uncolored diagrams in this case 
is equal to the number of matchings on b letters, and each of the b/2 edges of the diagram 
can be colored in one of a ways. Hence we have 

6< a h l 2 

ma)) = 

3.2.2. a even, b even. If a and b are even, we can allow an even number of isolated nodes. 
If our diagram contains i isolated nodes, the number of uncolored diagrams is equal to the 
number of permutations on b with shape {l l 2 b ~ l ), and each of the {b — i)/2 edges of the 
diagram can be colored in one of a ways. Hence we have 

b _ A b\a^/ 2 

i=0,i even \ 2 ' 

3.2.3. a even, b odd. Similarly, if b is odd we can allow an odd number of isolated nodes, 
and so 

b\a^/ 2 



=l,i odd v 2 



3.3. The General Case A = {l bl 2 b2 ■ ■ -t bt ). 

Lemma 4. Let g be a permutation and and let a be a matching such that ago — g. If (a (3) 
appears in a, then a, (5 € {1,2, ... ,n} appear in cycles of the same length in g. 

Proof. Suppose a appears in the cycle (aai a 2 ■ ■ • at&) and (3 appears in the cycle (/5 /?i /3 2 • • • A) 
Assume without loss of generality that k < I. We have 

aii = g(a) = crga(a) = ag(/3) = cx(/3i) 

and so (oti Pi) is a transposition appearing in a. Similarly, 

"2 = g(ai) = crga(ai) = ag(/3i) = a((3 2 ) 

and so (01% P2) is also a transposition appearing in a. Continuing in this way, we find that 
o must contain the transpositions 

(a 13), (ati (3i), (a 2 (3 2 ), ■■■,(a k 13k) 

Suppose for contradiction that k < I. Then there must exist some 7 G {1,2, ... ,n} such 
that 7 ^ {a, aii, ■ ■ ■ , «fc} and (7 (3i) is a transposition appearing in a. Then 

13 = g((3i) = crgo-((3i) = ag(>y) 
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By hypothesis (3 = a (a), so a = g(j)- But a = g(aik), a contradiction, since 7 7^ Hence 
k = l. □ 

Thus, if g has hi cycles of length i, we need only count the number of matchings which 
commute with the product of these bi cycles, which we can calculate using the formulas 
already determined. To find the total number of matchings which commute with g, we 
simply take the product. That is, 

N(X) = iV((l bl ))iV((2 62 )) • • • N((t bt )) 

As an example, consider the permutation 

g = (1 14 7 11)(5 16 12)(6 20 9) (8 17 13)(10 2 15 4)(18 3 19) 

We note that g has two cycles of length four, and by Lemma [], we have 

iV((4 2 )) = 4 + 1 = 5 

In addition, g has four cycles of length three, and 

4'3 4 / 2 

So the number of matchings that will commute with g is 5 ■ 27 = 135. 

3.4. Data. Table 1 displays dimension of V 2m (V(n)) Klyn ^ for several values of r and m 
generated by the formulas of Theorems || and |2|. 



r\m 


1 


2 


3 


4 


5 


6 


1 


1 


1 


1 


1 


1 


1 


2 


1 


2 


3 


5 


7 


11 


3 


1 


5 


16 


86 


448 


3580 


4 


1 


14 


132 


4154 


234004 


24791668 


5 


1 


41 


1439 


343101 


208796298 


253588562985 


6 


1 


122 


18373 


33884745 


196350215004 


2634094790963313 


7 


1 


365 


254766 


3505881766 


185471778824111 


27380169200102651288 


8 


1 


1094 


3680582 


366831842914 


175264150734326927 


284615877731708760168866 




Table 1: 


dimV 2m (V(n 


))A'(n) where K ( n j _ 


nLi OK), v» = riLi c n 


, and rii > 2m for 1 < j < r 



4. The Invariants and a Graph Interpretation 

We now seek a description of the invariants themselves. We begin by setting up notation 
for an arbitrary tensor in V^n). Let e» G C n denote the vector with a 1 in row i and 
elsewhere. Note that (ei, . . . , e n ) is an ordered basis for C n . An arbitrary tensor in V(n) 
is of the form 

where Zj ranges from 1 to nj and x^...^ is a complex scalar. 
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Recall that I\ m denotes the set of r-tuples of matchings on 2m letters. The group S^m 
acts on this set by simultaneous conjugation. We choose a representative (n, r 2 , . . . , r r ) from 
each orbit under this action, which we denote by 



Fix an ordering on the cycles of each Tj, and let j\ denote the cycle containing k in Tj. For 
instance, if r 2 = (1 3)(2 4), then j| = 2, since 4 appears in the second cycle of r 2 . The 
invariant associated to [ti, r 2 , . . . , r r ] can now written as the sum of terms of the form 



where each a\ ranges from 1 to n&. 

A /c-regular graph on n vertices is a graph in which each of the n vertices has degree k\ 
that is, each vertex is met by exactly k edges. An edge coloring of a graph is a coloring of 
these edges so that no two adjacent edges are the same color (this is also referred to as a 
1-factorization of the graph). In this case of a k- regular graph, an edge coloring implies that 
each vertex is met by k edges, which are each a distinct color. 

We now present a bijection between the orbits of the action of S2 m on 7£ m and the iso- 
morphism classes of edge colored r-regular graphs with 2m vertices. To construct the graph 
associated to [r 1; r 2 , . . . , r r ], number the vertices of the graph from 1 to 2m. An edge of color 
i is drawn between vertices j and k if Tj contains the cycle (j k). Repeating this process for 
all r matchings, we obtain an undirected graph with 2m vertices and r colors. 

To illustrate, we consider the case r = 3, m = 2 where we have chosen the 3-tuple of 
matchings 



To construct the graph associated to [ti, r 2 , r 3 ], we encode t±, r 2 , T3 with the colors black, red, 
and blue, respectively. Forgetting the labels of the vertices leaves us with a representative 
of the isomorphism class containing the graph. 



[ti, r 2 , . . . , r r ] = Win, r 2 , . . . , T r )(J I cr G S 2m } 




(Ti, t 2) t 3 ) = ((1 2)(3 4), (1 3) (2 4), (1 4) (2 3)) 



The invariant associated to [ti,t 2 ,T3] is then 




a[ r \a^ a< 2 \4 2 ) a^af 



Graph Q associated 



Representative of 
isomorphism class of Q 





Note that any graph isomorphic to this graph will correspond to a 3-tuple of matchings that 
belongs to the same orbit as our original choice. 
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A list of representatives of all isomorphism classes of 3-regular graphs on four vertices is 
shown in Table 2, along with the corresponding invariant. We have chosen a representative 
(ti,t 2 ,t 3 ) of each orbit so that t\ = (1 2) (3 4). 

Finally, we show a way of encoding these invariants with forests of phylogenetic trees. 
As explored in the introduction, Diaconis and Holmes | DH98 | provide a bijection between 
matchings on 2m letters and phylogenetic trees with m + 1 labelled leaves. When given an 
r-tuple of matchings, we can define a forest of r trees with roots labelled 1,2, ... ,r. Note 
that we consider two forests to be equivalent if the individual trees within the forest are 
equivalent in the traditional way. As an example, we again consider the 3-tuple 



(ri,r 2 ,r 3 ) = ((1 2) (3 4), (1 3)(2 4), (1 4) (2 3)) 



which can be interpreted as the forest 




after we have labelled all vertices following the method outlined in section 1. 

The invariant associated to a forest with r trees can again be written as a sum of terms 
of the form 

X (1) (2) (r)X (1) (2) (r) ■ • • X (1) (2) (r) 

a , a ., —a , a „ a , ■■■a\ a V™ a V™ —a \L. 

.■1 .'1 ^-1 ,;2 ,-2 a2 ,-2m ^2m „-2m 

1 2 r 1 2 r -"1 ^2 ■' r 

where the index (i) of the subscript refers to the root of an individual tree in the forest, and 
it = j\ if s an d t are siblings in the ith tree. To build the invariant associated with our 
particular forest, we begin by examining the tree labelled (1). This tree has two pairs of 
siblings, (1,2) and (3,4), and so our invariant begins as 

X (i)__ X (i)__ X (i)__ X (i)__ 

Next, we see that the second tree has sibling pairs (1,3) and (2,4), so we continue building 
the invariant: 



X (1) (2)_X (1) (2)-X (1) (2)_X (1) (2). 



Finally, the last placeholders of the invariant are filled by observing that the third tree has 
sibling pairs (2,3) and (1,4). The result is identical to the invariant associated to the same 
tuple of matchings determined earlier. 



X (1) (2) (3)X (1) (2) (3)X (1) (2) (3)X (1) (2) (3) 

14 





V^ n l V^ n 2 V^ n 3 

2^(1) (1)2^(2) (2)2^(3) (3) X „m n {2) n {3) X n {l) n {2) n (3) X „(1) „(2) „(3)X (1) (2) (3) 
,aj ,a;, a\ ' ,a 2 «i o- l a 1 a 1 a l a l a 2 a 2 a 2 a 2 a 2 a 2 








V^ ra l V^ n 2 V^ n 3 

2^(1) (1)2^(2) (2)2^(3) (3) £„(1)„(2)„(3)£ „(2) (3)X (1) (2) (3)X (1) (2) (3) 
a x ,a 2 a\ ,a 2 a\ ,a 2 % «i «i % «i «2 a 2 a 2 a 2 a 2 a 2 a i 








v-^ni v^ n 2 V^"3 

lu (1) (1)2^(2) (2)2^(3) (3) ^„(1).(2)„(3)X (1) (2) (3)X (1) (2) (3)X (1) (2) (3) 
jOj ,aj Oj jOj a l a l a l a l a 2 a l a 2 a 2 a 2 a 2 a l a 2 





\—yni V^ n 2 V^™3 

2^(1) (1)2^(2) (2)2^(3) &X(i) (2) (Z)X(l) <2) (3)X(l) (2) (i)X(i) (2) (i) 
a\ ,a\ a\ ' ,a 2 a\ ' ,a> 2 a 1 a± a± a 1 a 2 a 2 a 2 a 2 a 2 a 2 a t a x 






> 


V^ ra l v^ n 2 V^"3 

2^(1) (1)2^(2) (2)2^(3) (3) X „( 1 ) f ,( 2 )„( 3 ) ;r „( 1 )„( 2 )„( 3 ) ;r „( 1 )„( 2 )„( 3 ) ;Z; „( 1 )„( 2 )„( 3 ) 
*-^a\ ,a 2 ' ,a\ 'a\ ,a 2 a 1 a 1 a 1 a 1 a 2 a 2 a 2 a 1 a 2 a 2 a 2 a 1 



Table 2: Representative of graph isomorphism class and corresponding invariant in the case 
r = 3, m = 2 
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